Precision Measurement of the Neutron Spin Asymmetries and Spin-dependent Structure Functions 

in the Valence Quark Region 
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We report on measurements of the neutron spin asymmetries A™ 2 an d polarized structure functions g™ 2 at 
three kinematics in the deep inelastic region, with x = 0.33, 0.47 and 0.60 and Q 2 = 2.7, 3.5 and 4.8 (GeV/c) 2 , 
respectively. These measurements were performed using a 5.7 GeV longitudinally-polarized electron beam and 
a polarized 3 He target. The results for A" and g" at x = 0.33 are consistent with previous world data and, at 
the two higher x points, have improved the precision of the world data by about an order of magnitude. The 
new A" data show a zero crossing around x = 0.47 and the value at x = 0.60 is significantly positive. These 
results agree with a next-to-leading order QCD analysis of previous world data. The trend of data at high x 
agrees with constituent quark model predictions but disagrees with that from leading-order perturbative QCD 
(pQCD) assuming hadron helicity conservation. Results for A 2 and g 2 have a precision comparable to the best 
world data in this kinematic region. Combined with previous world data, the moment d 2 was evaluated and 
the new result has improved the precision of this quantity by about a factor of two. When combined with the 
world proton data, polarized quark distribution functions were extracted from the new g" / F" values based on 
the quark parton model. While results for Au/u agree well with predictions from various models, results for 
Ad/d disagree with the leading-order pQCD prediction when hadron helicity conservation is imposed. 



PACS numbers: 13.60.Hb,24.85.+p,25.30.-c 
I. INTRODUCTION 

Interest in the spin structure of the nucleon became promi- 
nent in the 1980's when experiments at CERN |jl|] and 



SLAC 01 on the integral of the proton polarized structure 
function g\ showed that the total spin carried by quarks was 
very small, w (12±17)% d. This was in contrast to the sim- 
ple relativistic valence quark model prediction j3| in which 
the spin of the valence quarks carries approximately 75% of 
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the proton spin and the remaining 25% comes from their or- 
bital angular momentum. Because the quark model is very 
successful in describing static properties of hadrons, the fact 
that the quark spins account for only a small part of the nu- 
cleon spin was a big surprise and generated very productive 
experimental and theoretical activities to the present. Current 
understanding @] of the nucleon spin is that the total spin is 
distributed among valence quarks, qq sea quarks, their orbital 
angular momenta, and gluons. This is called the nucleon spin 
sum rule: 



S 



N 



1 



(1) 



where S% is the nucleon spin, 5? and L\ represent respec- 
tively the quark spin and orbital angular momentum (OAM), 
and J| is the total angular momentum of the gluons. Only 
about (20 — 30)% of the nucleon spin is carried by the spin of 
the quarks. To further study the nucleon spin, one thus needs 
to know more precisely how it decomposes into the three com- 
ponents and to measure their dependence on x. Here x is the 
Bjorken scaling variable, which in the quark-parton model 
can be interpreted as the fraction of the nucleon momentum 
carried by the quark. For a fixed target experiment one has 
x = Q 2 1 (2Mv), with M the nucleon mass, Q 2 the four mo- 
mentum transfer squared and v the energy transfer from the 
incident electron to the target. However, due to experimen- 
tal limitations, precision data have been collected so far only 
in the low and moderate x regions. In these regions, one is 
sensitive to contributions from a large amount of qq sea and 
gluons and the nucleon is hard to model. Moreover, at large 
distances corresponding to the size of a nucleon, the theory 
of the strong interaction - Quantum Chromodynamics (QCD) 
- is highly non-perturbative, which makes the investigation 
of the roles of quark orbital angular momentum (OAM) and 
gluons in the nucleon spin structure difficult. 

Our focus here is the first precise neutron spin structure data 
in the large x region x > 0.4. For these kinematics, the va- 
lence quarks dominate and the ratios of structure functions 
can be estimated based on our knowledge of the interactions 
between quarks. More specifically, the virtual photon asym- 
metry A\, defined as 



Ai{x, 
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01/2 + 03/2 



(the definitions of 0-1/2,3/2 are given in Appendix lA"li. which 
at large Q 2 is approximately the ratio of the polarized and 
the unpolarized structure functions g\/F\, is expected to 
approach unity as x — * 1 in perturbative QCD (pQCD). This 
is a dramatic prediction, not only because this is the only 
kinematic region where one can give an absolute prediction 
for the structure functions based on pQCD, but also because 
all previous data on the neutron asymmetry A r { in the region 
x > 0.4 have large uncertainties and are consistent with 
A™ ^ 0. Furthermore, because both qq sea and gluon 
contributions are small in this region, it is a relatively clean 
region to test the valence quark model and to study the role 
of valence quarks and their OAM contribution to the nucleon 



spin. 



Deep inelastic scattering (DIS) has served as one of the ma- 
jor experimental tools to study the quark and gluon structure 
of the nucleon. The formalism of unpolarized and polarized 
DIS is summarized in Appendix [X] Within the quark par- 
ton model (QPM), the nucleon is viewed as a collection of 
non-interacting, point-like constituents, one of which carries 
a fraction x of the nucleon's longitudinal momentum and ab- 
sorbs the virtual photon @|. The nucleon cross section is then 
the incoherent sum of the cross sections for elastic scattering 
from individual charged point-like partons. Therefore the un- 
polarized and the polarized structure functions Fx and gi can 
be related to the spin-averaged and spin-dependent quark dis- 
tributions as 1 6 1 



and 
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(x,Q 2 )=i^ e 2A g! :(x,Q 2 ), 



(2) 



(3) 



where qi(x, Q 2 ) = qj (x, Q 2 ) + q\(x, Q 2 ) is the unpolarized 
parton distribution function (PDF) of the i th quark, defined 
as the probability that the i th quark inside a nucleon carries 
a fraction x of the nucleon's momentum, when probed with a 
resolution determined by Q 2 . The polarized PDF is defined 
as Aq t (x,Q 2 ) = qj(x,Q 2 ) - q\{x,Q 2 \ where qj(x,Q 2 ) 
(q\(x, Q 2 )) is the probability to find the spin of the i th quark 
aligned parallel (anti-parallel) to the nucleon spin. 

The polarized structure function 92 (x, Q 2 ) does not have a 
simple interpretation within the QPM |6]. However, it can be 
separated into leading twist and higher twist terms using the 
operator expansion method 01 : 



(4) 



Here g\ vw {x : Q 2 ) is the leading twist (twist-2) contribu- 
tion and can be calculated using the twist-2 component of 
gi(x, Q 2 ) and the Wandzura-Wilczek relation 0] as 



f 1 9i(y, 



%■ (5) 



g™ w {x,Q 2 ) = - gi (x,Q 2 ) + 

j x y 

The higher-twist contribution to 92 is given by g 2 . When ne- 
glecting quark mass effects, the higher-twist term represents 
interactions beyond the QPM, e.g., quark-gluon and quark- 
quark correlations The moment of g 2 can be related to the 
matrix element d 2 11011 : 



da; 



la 



3g 2 (x 1 Q 2 ) + 2g 1 (x,Q 2 ) 



3 / dx x 2 g 2 (x, Q 2 ) 
Jo 



(6) 



Hence d 2 measures the deviations of g 2 from g\ vw . The value 
of d 2 can be obtained from measurements of g\ and g 2 and 
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can be compared with predictions from Lattice QCD fiUl . 
bag models Il2ll . QCD sum rules lU3ll and chiral soliton mod- 
els E~ 

In this paper we first describe available predictions for A\ 
at large x. The experimental apparatus and the data analysis 
procedure will be described in Section lTTlllIVI andlVl In Sec- 
tion^] we present results for the asymmetries and polarized 
structure functions for both 3 He and the neutron, a new exper- 
imental fit for gf/FJ 1 and a result for the matrix element d%. 
Combined with the world proton and deuteron data, polarized 
quark distribution functions were extracted from our / F\ 
results. We conclude the paper by summarizing the results for 
A™ and Ad/d and speculating on the importance of the role 
of quark OAM on the nucleon spin in the kinematic region 
explored. Some of the results presented here were published 
previously fl5ll : the present publication gives full details on 
the experiment and all of the neutron spin structure results for 
completeness. 



II. PREDICTIONS FOR A" AT LARGE x 

From Section Hi Al to IITfI we present predictions of A" at 
large x. Data on A" from previous experiments did not have 
the precision to distinguish among different predictions, as 
will be shown in SectionlHGl 



proton, s n (x) = s p (x) = s(x). In the following, all PDF's 
are for the proton, unless specified by a superscript 'n'. 

In the case of DIS, exact SU(6) symmetry implies the same 
shape for the valence quark distributions, i.e. u(x) = 2d(x). 
Using Eq. (0 and iA4\ . and assuming that R(x, Q 2 ) is the 
same for the neutron and the proton, one can write the ratio of 
neutron and proton F2 structure functions as 



R np 



n 



u(x) + 4d(x) 



4u(x) + d{x) 



(10) 



Applying u(x) = 2d(x) gives 



R np = 2/3. (11) 

However, data on the R np ratio from SLAC O, CERN (H 
Eel I2II1 and Fermilab I22I1 disagree with this SU(6) predic- 
tion. The data show that R np (x) is a straight line starting with 
R np \x^o ~ 1 and dropping to below 1/2 as x — > 1. In ad- 
dition, A\(x) is small at low x Q IS Q ■ The fact that 
R np \ x ^o « 1 may be explained by the presence of a domi- 
nant amount of sea quarks in the low x region and the fact that 
A\ | x ^o ~ could be because these sea quarks are not highly 
polarized. At large x, however, there are few sea quarks and 
the deviation from SU(6) prediction indicates a problem with 
the wavefunction described by Eq. Q. In fact, SU(6) sym- 
metry is known to be broken 1 26] and the details of possi- 
ble SU(6)-breaking mechanisms is an important open issue in 
hadronic physics. 



A. SU(6) Symmetric Non-Relativistic Constituent Quark 
Model 

In the simplest non-relativistic constituent quark 
model (CQM) lUal . the nucleon is made of three con- 
stituent quarks and the nucleon spin is fully carried by the 
quark spin. Assuming SU(6) symmetry, the wavefunction of a 
neutron polarized in the +z direction then has the form fl7ll : 

\n T) = I^Wooo) + -4= \d\du) lw ) (7) 
y 2 v 18 

1 1 

-- \d l (du)m) - - |w T (<id)iio) + — \u l (dd)m) , 

where the three subscripts are the total isospin, total spin 5 
and the spin projection S z along the +z direction for the 'di- 
quark' state. For the case of a proton one needs to exchange 
the u and d quarks in Eq. 0. In the limit where SU(6) sym- 
metry is exact, both diquark spin states with 5 = 1 and 5 = 
contribute equally to the observables of interest, leading to the 
predictions 

A\ = 5/9 and A" = ; (8) 
Au/u -> 2/3 and Ad/d -> -1/3. (9) 

We define u{x) = u p (x), d(x) = d p (x) and s(x) = s p (x) 
as parton distribution functions (PDF) for the proton. For a 
neutron one has u n (x) = d p (x) = d(x), d n (x) = u p (x) = 
u{x) based on isospin symmetry. The strange quark distribu- 
tion for the neutron is assumed to be the same as that of the 



B. SU(6) Breaking and Hyperfine Perturbed Relativistic CQM 

A possible explanation for the SU(6) symmetry breaking 
is the one-gluon exchange interaction which dominates the 
quark-quark interaction at short-distances. This interaction 
was used to explain the behavior of R np near x — * 1 and 
the w 300-MeV mass shift between the nucleon and the 
A(1232) l2(ill . Later this was described by an interaction term 
proportional to Si ■ Sj S 3 (fij), with 5; the spin of the i th 
quark, hence is also called the hyperfine interaction, or chro- 
momagnetic interaction among the quarks l27il . The effect of 
this perturbation on the wavefunction is to lower the energy 
of the 5 = diquark state, causing the first term of Eq. @, 
\d t (ud)ooo) n , to become more stable and to dominate the 
high energy tail of the quark momentum distribution that is 
probed as x — > 1. Since the struck quark in this term has its 
spin parallel to that of the nucleon, the dominance of this term 
asi-»l implies (Ad/d) n -> 1 and (Au/u) n -> -1/3 for 
the neutron, while for the proton one has 

Au/u -»■ 1 and Ad/d-> -1/3 asx -> 1 . (12) 

One also obtains 

R np -> 1/4 asx -> 1 , (13) 

which could explain the deviation of R np (x) data from the 
SU(6) prediction. Based on the same mechanism, one can 
make the following predictions: 

A\—>1 and A™ — > 1 as x -> 1 . (14) 
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The hyperfine interaction is often used to break SU(6) sym- 
metry in the relativistic CQM (RCQM). In this model, the 
constituent quarks have non-zero OAM which carries w 25% 
of the nucleon spin 01 . The use of RCQM to predict the large 
x behavior of the nucleon structure functions can be justified 
by the valence quark dominance, i.e., in the large x region 
almost all quantum numbers, momentum and the spin of the 
nucleon are carried by the three valence quarks, which can 
therefore be identified as constituent quarks. Predictions of 
A™ and A p in the large x region using the hyperfine-perturbed 
RCQM have been achieved lEil . 

C. Perturbative QCD and Hadron Helicity Conservation 

In the early 1970's, in one of the first applications of per- 
turbative QCD (pQCD), it was noted that as a; — ► 1, the scat- 
tering is from a high-energy quark and thus the process can 
be treated perturbatively 12911 . Furthermore, when the quark 
OAM is assumed to be zero, the conservation of angular mo- 
mentum requires that a quark carrying nearly all the momen- 
tum of the nucleon (i.e. x — ► 1) must have the same helicity as 
the nucleon. This mechanism is called hadron helicity conser- 
vation (HHC), and is referred to as the leading-order pQCD 
in this paper. In this picture, quark-gluon interactions cause 
only the S = 1, S z = 1 diquark spin projection component 
rather than the full 5 = 1 diquark system to be suppressed as 
x — > 1, which gives 

Au/u — > 1 and Ad/d — > 1 as x — > 1 ; (15) 

R np -, A\ -»• 1 and A\ -> 1 as x -> 1 . (16) 

This is one of the few places where pQCD can make an abso- 
lute prediction for the x-dependence of the structure functions 
or their ratios. However, how low in x and Q 2 this picture 
works is uncertain. HHC has been used as a constraint in a 
model to fit data on the first moment of the proton g\, giv- 
ing the BBS parameterization iBOll . The Q 2 evolution was not 
included in this calculation. Later in the LSS(BBS) parame- 
terization l3lll . both proton and neutron A\ data were fitted di- 
rectly and the Q 2 evolution was carefully treated. Predictions 
for A™ using both BBS and LSS(BBS) parameterizations have 
been made, as shown in Fig.^ an d|2]in Section HTGl 

HHC is based on the assumption that the quark OAM 
is zero. Recent experimental data on the tensor polariza- 
tion in elastic e— 2 H scattering l32ll . neutral pion photo- 
production l33ll and the proton electro-magnetic form fac- 
tors O Q disagree with the HHC predictions It has 
been suggested that effects beyond leading-order pQCD, such 
as quark OAM 1371 l38i l39l 14011 . might play an important role 
in processes involving quark spin flips. 

D. Predictions from Next-to-Leading Order QCD Fits 

In a next-to-leading order (NLO) QCD analysis of the 
world data l4lll . parameterizations of the polarized and un- 
polarized PDFs were performed without the HHC constraint. 



Predictions of g p / Ff and / were made using these pa- 
rameterizations, as shown in Fig.^andOlin Section HTGl 

In a statistical approach, the nucleon is viewed as a gas of 
massless partons (quarks, antiquarks and gluons) in equilib- 
rium at a given temperature in a finite volume, and the parton 
distributions are parameterized using either Fermi-Dirac or 
Bose-Einstein distributions. Based on this statistical picture 
of the nucleon, a global NLO QCD analysis of unpolarized 
and polarized DIS data was performed 14211 . In this calcula- 
tion Au/u w 0.75, Ad/d w -0.5 and A\' n < 1 at x -> 1. 

E. Predictions from Chiral Soliton and Instanton Models 

While pQCD works well in high-energy hadronic physics, 
theories suitable for hadronic phenomena in the non- 
perturbative regime are much more difficult to construct. Pos- 
sible approaches in this regime are quark models, chiral effec- 
tive theories and the lattice QCD method. Predictions for A™' p 
have been made using chiral soliton models I43I I44I and the 
results of Ref. Q give A™ < 0. The prediction that A\ < 
has also been made in the instanton model 14511 . 



F. Other Predictions 

Based on quark-hadron duality l46Tl . one can obtain the 
structure functions and their ratios in the large x region by 
summing over matrix elements for nucleon resonance tran- 
sitions. To incorporate SU(6) breaking, different mechanisms 
consistent with duality were assumed and data on the structure 
function ratio R np were used to fit the SU(6) mixing param- 
eters. In this picture, A™ :P — > 1 as x — ► 1 is a direct result. 
Duality predictions for A™' p using different SU(6) breaking 
mechanisms were performed in Ref. l47Tl . There also exist 
predictions from bag models 0, as shown in Fig.^ an dElin 
the next section. 



G. Previous Measurements of A" 

A summary of previous A™ measurements is given in 



TABLE I: Previous measurements of A". 



Experiment 


beam 


target 


X 


Q 2 

(GeV/c) 2 


E142 [5f] 


19.42, 22.66, 


3 He 


0.03-0.6 


2 


25.51 GeV; e" 








E154[52J 


48.3 GeV; e" 


3 He 


0.014-0.7 


1-17 


HERMES [50J 


27.5 GeV; e+ 


3 He 


0.023-0.6 


1-15 


E143 [251 


9.7, 16.2, 


NH 3 , ND 3 


0.024-0.75 


0.5-10 




29.1 GeV; e~ 








E155 [53J 


48.35 GeV; e~ 


NH 3 ,LiD 3 


0.014-0.9 


1-40 


SMC [491 


190 GeV; /x" 


C 4 Hi„0 


0.003-0.7 


1-60 




C4D10O 
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0.2 0.4 0.6 0.8 1 

X 

FIG. 1: Previous data on A? |H El |5(J El |H |H and vari- 
ous theoretical predictions: A" from SU(6) symmetry (solid line at 
zero) [n|], hyperfine-perturbed RCQM (shaded bandVjm BBS pa- 
rameterization at Q 2 = 4 (GeV/c) 2 (higher solid) 1301. LSS(BBS) 
parameterization at Q 2 = 4 (GeV/c) 2 (dashed) 13111 . statistical model 
at Q 2 — A (GeV/c) 2 (long-dashed) 1421 . quark-hadron duality us- 
ing two different SU(6) breaking mechanisms (dash-dot-dotted and 
dash-dot-dot-dotted) [47], and non-meson cloudy bag model (dash- 
dotted) EH; g?/F? from LSS2001 parameterization at Q 2 = 
5 (GeV/c) 2 (lower solid) l4lll and from chiral soliton models l43Tl at 
Q 2 = 3 (GeV/c) 2 (long dash-dotted) and Q at g 2 = 4.8 (GeV/c) 2 
(dotted). 

Table U The data on A™ and Aj 1 are plotted in Fig. ^ and|2| 
along with theoretical calculations described in previous sec- 
tions. Since the Q 2 -dependence of A\ is small and g\jF\ k, 
A\ in DIS, data for gf /F[ l and g\ j Ff are also shown and all 
data are plotted without evolving in Q 2 . As becomes obvi- 
ous in Fig.[U the precision of previous A™ data at x > 0.4 
from SMC [49], HERMES H3] and SLAC HUB HI is not 
sufficient to distinguish among different predictions. 



III. THE EXPERIMENT 

We report on an experiment l55ll carried out at in the Hall A 
of Thomas Jefferson National Accelerator Facility (Jefferson 
Lab, or JLab). The goal of this experiment was to provide 
precise data on A™ in the large x region. We have measured 
the inclusive deep inelastic scattering of longitudinally polar- 
ized electrons off a polarized 3 He target, with the latter being 
used as an effective polarized neutron target. The scattered 
electrons were detected by the two standard High Resolution 
Spectrometers (HRS). The two HRS were configured at the 
same scattering angles and momentum settings to double the 
statistics. Data were collected at three x points as shown in 
Table UTI Both longitudinal and transverse electron asymme- 




0.2 0.4 0.6 0.8 1 

X 



FIG. 2: World data on A\ O ^ |3 S IH and predictions for 
gl/Fl at Q 2 = 5 (GeV/c) 2 from the E155 experimental fit (long 
dash-dot-dotted) l53ll and a new fit as described in Section fVI Bl flong 
dash-dot-dot-dotted). The solid curve corresponds to the prediction 
for g?/F? from LSS(2001) parameterization at Q 2 = 5 (GeV/c) 2 . 
Other curves are the same as in Fig.Qexcept that there is no predic- 
tion for the proton from BBS and LSS(BBS) parameterizations. 



tries were measured, from which A\, A^, gi/Fi and gijFx 
were extracted using Eq. JA22tjA"25l . 

TABLE II: Kinematics of the experiment. The beam energy was 
E = 5.734 GeV. E 1 and 6 are the nominal momentum and angle 
of the scattered electrons, (x), (Q 2 ) and (W 2 ) are values averaged 
over the spectrometer acceptance. 



(x) 


0.327 


0.466 


0.601 


E' 

e 

(Q 2 ) (GeV/c) 2 
{W 2 } (GeV) 2 


1.32 
35° 
2.709 
6.462 


1.72 
35° 
3.516 
4.908 


1.455 
45° 
4.833 
4.090 



A. Polarized 3 He as an Effective Polarized Neutron 

As shown in Fig.^ previous data on A™ did not have suf- 
ficient precision in the large x region. This is mainly due to 
two experimental limitations. Firstly, high polarization and 
luminosity required for precision measurements in the large x 
region were not available previously. Secondly, there exists no 
free dense neutron target suitable for a scattering experiment, 
mainly because of the neutron's short lifetime (« 886 sec). 
Therefore polarized nuclear targets such as 2 H or 3 Hc are 
commonly used as effective polarized neutron targets. Con- 
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sequently, nuclear corrections need to be applied to extract 
neutron results from nuclear data. For a polarized deuteron, 




S state D state S' state 

(89.93%) (8.75%) (1.26%) 



FIG. 3: An illustration of 5 He wavefunction. The S, S' and D state 
contributions are from calculations using the AVI 8 two-nucleon in- 
teraction and the Tucson-Melbourne three-nucleon force, as given in 
Ref. IH. 

approximately half of the deuteron spin comes from the pro- 
ton and the other half comes from the neutron. Therefore the 
neutron results extracted from the deuteron data have a signifi- 
cant uncertainty coming from the error in the proton data. The 
advantage of using 3 Hc is that the two protons' spins cancel 
in the dominant S state of the 3 He wavefunction, thus the spin 
of the 3 He comes mainly (> 87%) from the neutron 1561 15711 . 
as illustrated in Fig. || As a result, there is less model de- 
pendence in the procedure of extracting the spin-dependent 
observables of the neutron from 3 He data. At large x, the ad- 
vantage of using a polarized 3 He target is more prominent in 
the case of A™. In this region almost all calculations show 
that A\ is much smaller than A\, therefore the A" results ex- 
tracted from nuclear data are more sensitive to the uncertainty 
in the proton data and the nuclear model being used. 

In the large x region, the cross sections are small because 
the parton densities drop dramatically as x increases. In ad- 
dition, the Mott cross section, given by Eq. IA3I is small at 
large Q 2 . To achieve a good statistical precision, high lumi- 
nosity is required. Among all laboratories which are equipped 
with a polarized 3 He target and are able to perform a mea- 
surement of the neutron spin structure, the polarized electron 
beam at JLab, combined with the polarized 3 He target in Hall 
A, provides the highest polarized luminosity in the world lEH . 
Hence it is the best place to study the large x behavior of the 
neutron spin structure. 



B. The Accelerator and the Polarized Electron Source 

JLab operates a continuous-wave electron accelerator that 
recirculates the beam up to five times through two super- 
conducting linear accelerators. Polarized electrons are ex- 
tracted from a strained GaAs photocathode |59] illuminated 
by circularly polarized light, providing a polarized beam of 
(70 — 80)% polarization and « 200/iA maximum current to 
experimental halls A, B and C. The maximum beam energy 
available at JLab so far is 5.7 GeV, which was also the beam 
energy used during this experiment. 



C. Hall A Overview 

The basic layout of Hall A during this experiment is shown 
in Fig. 0] The major instrumentation l60ll includes beamline 
equipment, the target and two HRSs. The beamline starts af- 




FIG. 4: (Color online) Top-view of the experimental hall A (not to 
scale). 

ter the arc section of the accelerator where the beam is bent 
into the hall, and ends at the beam dump. The arc section 
can be used for beam energy measurement, as will be de- 
scribed in Section ITlI Dl After the arc section, the beamline 
is equipped with a Compton polarimeter, two Beam Current 
Monitors (BCM) and an Unser monitor for absolute beam cur- 
rent measurement, a fast raster, the eP device for beam en- 
ergy measurement, a M0ller polarimeter and two Beam Posi- 
tion Monitors (BPM). These beamline elements, together with 
spectrometers and the target, will be described in detail in the 
following sections. 

D. Beam Energy Measurement 

The energy of the beam was measured absolutely by two 
independent methods - ARC and eP l60l l6lll . Both methods 
can provide a precision of 8Eb eam / Eb eam ps 2 x 10~ 4 . For 
the ARC method |6(J|62|], the deflection of the beam in the arc 
section of the beamline is used to determine the beam energy. 
In the eP measurement ^3^3l the beam energy is determined 
by the measurement of the scattered electron angle 9 e and the 
recoil proton angle 6 p in 1 H(e, e'p) elastic scattering. 

E. Beam Polarization Measurement 

Two methods were used during this experiment to measure 
the electron beam polarization. The M0ller polarimeter ll60ll 
measures M0ller scattering of the polarized electron beam off 
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polarized atomic electrons in a magnetized foil. The cross 
section of this process depends on the beam and target po- 
larizations. The polarized electron target used by the M0ller 
polarimeter was a ferromagnetic foil, with its polarization de- 
termined from foil magnetization measurements. The M0ller 
measurement is invasive and typically takes an hour, provid- 
ing a statistical accuracy of about 0.2%. The systematic error 
comes mainly from the error in the foil target polarization. 
An additional systematic error is due to the fact that the beam 
current used during a M0ller measurement (« 0.5/xA) is lower 
than that used during the experiment. The total relative sys- 
tematic error was w 3.0% during this experiment. 

During a Compton polarimeter l60l l64ll measurement, the 
electron beam is scattered off a circularly polarized photon 
beam and the counting rate asymmetry of the Compton scat- 
tered electrons or photons between opposite beam helicities 
is measured. The Compton polarimeter measures the beam 
polarization concurrently with the experiment running in the 
hall. 

The Compton polarimeter consists of a magnetic chicane 
which deflects the electron beam away from the scattered pho- 
tons, a photon source, an electromagnetic calorimeter and an 
electron detector. The photon source was a 200 mW laser am- 
plified by a resonant Fabry-Perot cavity. During this experi- 
ment the maximum gain of the cavity reached G max = 7500, 
leading to a laser power of 1500 W inside the cavity. The 
circular polarization of the laser beam was > 99% for both 
right and left photon helicity states. The asymmetry measured 
in Compton scattering at JLab with a 1.165 eV photon beam 
and the 5.7 GeV electron beam used by this experiment had a 
mean value of ~ 2.2% and a maximum of 9.7%. For a 12 /iA 
beam current, one hour was needed to reach a relative statis- 
tical accuracy of (AP&) s t a t I Pb ~ 1%. The total systematic 
error was (APf,) sys / Pj, w 1.6% during this experiment. 

The average beam polarization during this experiment was 
extracted from a combined analysis of 7 M0ller and 53 Comp- 
ton measurements. A value of (79.7 ± 2.4)% was used in the 
final DIS analysis. 



F. Beam Helicity 

The helicity state of electrons is regulated every 33 ms at the 
electron source. The time sequence of the electrons' helicity 
state is carried by helicity signals, which are sent to exper- 
imental halls and the data acquisition (DAQ) system. Since 
the status of the helicity signal (H+ or H- pulses) has either 
the same or the opposite sign as the real electron helicity, the 
absolute helicity state of the beam needs to be determined by 
other methods, as will be described later. 

There are two modes - toggle and pseudorandom - which 
can be used for the pulse sequence of the helicity signal. In 
the toggle mode, the helicity alternates every 33 ms. In the 
pseudorandom mode, the helicity alternates randomly at the 
beginning of each pulse pair, of which the two pulses must 
have opposite helicities in order to equalize the numbers of the 
H+ and H- pulses. The purpose of the pseudorandom mode is 
to minimize any possible time-dependent systematic errors. 



Fig. |5] shows the helicity signals and the helicity states of the 



DAQ + - + - + - + - + - + 



one pulse pair 



the helicity alternates randomly 
between pulse pairs 



33 ms 



DAQ 



FIG. 5: Helicity signal and the helicity status of DAQ in toggle (top) 
and pseudorandom (bottom) modes. 

DAQ system for the two regulation modes. 

There is a half-wave plate at the polarized source which 
can be inserted to reverse the helicity of the laser illuminating 
the photocathode hence reverse the helicity of electron beam. 
During the experiment this half-wave plate was inserted for 
half of the statistics to minimize possible systematic effects 
related to the beam helicity. 

The scheme described above was used to monitor the rela- 
tive changes of the helicity state. The absolute sign of the elec- 
trons' helicity states during each of the H+ and H- pulses were 
confirmed by measuring a well known asymmetry and com- 
paring the measured asymmetry with its prediction, as will be 
presented in Section lVBl and rVCl 



G. Beam Charge Measurement and Charge Asymmetry 
Feedback 



The beam current was measured by the BCM system lo- 
cated upstream of the target on the beamline. The BCM sig- 
nals were fed to scaler inputs and were inserted in the data 
stream. 

Possible beam charge asymmetry measured at Hall A can 
be caused by the timing asymmetry of the DAQ system, or 
by the timing and the beam intensity asymmetries at the po- 
larized electron source. The beam intensity asymmetry orig- 
inates from the intensity difference between different helicity 
states of the circularly polarized laser used to strike the pho- 
tocathode. Although the charge asymmetry can be corrected 
for to first order, there may exist unknown non-linear effects 
which can cause a systematic error in the measured asymme- 
try. Thus the beam charge asymmetry should be minimized. 
This was done by using a separate DAQ system initially devel- 
oped for the parity-violation experiments |65], called the par- 
ity DAQ. The parity DAQ used the measured charge asymme- 
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try in Hall A to control the orientation of a rotatable half-wave 
plate located before the photocathode at the source, such that 
intensities for each helicity state of the polarized laser used to 
strike the photocathode were adjusted accordingly. The parity 
DAQ was synchronized with the two HRS DAQ systems so 
that the charge asymmetry in the two different helicity states 
could be monitored for each run. The charge asymmetry was 
typically controlled to be below 2 x 10~ 4 during this experi- 
ment. 



H. Raster and Beam Position Monitor 

To protect the target cell from being damaged by the effect 
of beam-induced heating, the beam was rastered at the target. 
The raster consists a pair of horizontal and vertical air-core 
dipoles located upstream of the target on the beamline, which 
can produce either a rectangular or an elliptical pattern. We 
used a raster pattern distributed uniformly over a circular area 
with a radius of 2 mm. 

The position and the direction of the beam at the target 
were measured by two BPMs located upstream of the tar- 
get IrJoll . The beam position can be measured with a precision 
of 200 yum with respect to the Hall A coordinate system. The 
beam position and angle at the target were recorded for each 
event. 



I. High Resolution Spectrometers 

The Hall A High Resolution Spectrometer (HRS) systems 
were designed for detailed investigations of the structure of 
nuclei and nucleons. They provide high resolution in mo- 
mentum and in angle reconstruction of the reaction prod- 
uct as well as being able to be operated at high luminos- 
ity. For each spectrometer, the vertically bending design in- 
cludes two quadrupoles followed by a dipole magnet and a 
third quadrupole. All quadrupoles and the dipole are super- 
conducting. Both HRSs can provide a momentum resolution 
better than 2 x 10 -4 and a horizontal angular resolution bet- 
ter than 2 mrad with a design maximum central momentum 
of 4 GeV/c l60ll . By convention, the two spectrometers are 
identified as the left and the right spectrometers based on their 
position when viewed looking downstream. 

The basic layout of the left HRS is shown in Fig. [6] The 
detector package is located in a large steel and concrete 
detector hut following the last magnet. For this experiment 
the detector package included (1) two scintillator planes SI 
and S2 to provide a trigger to activate the DAQ electronics; 
(2) a set of two Vertical Drift Chambers (VDC) S for 
particle tracking; (3) a gas Cerenkov detector to provide 
particle identification (PID) information; and (4) a set of lead 
glass counters for additional PID. The layout of the right HRS 
is almost identical except a slight difference in the geometry 
of the gas Cerenkov detector and the lead glass counters. 
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FIG. 6: (Color online) Schematic layout of the left HRS and detector 
package (not to scale). 



J. Particle Identification 

For this experiment the largest background came from 
photo-produced pions. We refer to PID in this paper as the 
identification of electrons from pions. PID for each HRS was 
accomplished by a CO2 threshold gas Cerenkov detector and 
a double-layered lead glass shower detector. 

The two Cerenkov detectors, one on each HRS, were oper- 
ated with CO2 at atmospheric pressure. The refraction index 
of the CO2 gas was 1.00041, giving a threshold momentum of 
« 17 MeV/c for electrons and « 4.8 GeV/c for pions. The 
incident particles on each HRS were also identified by their 
energy deposits in the lead glass shower detector. 

Since Cerenkov detectors and lead glass shower detectors 
are based on different mechanisms and their PID efficiencies 
are not correlated Ifr/l . we extracted the PID efficiency of the 
lead glass counters by using electron events selected by the 
Cerenkov detector, and vice versa. Fig.0shows a spectrum of 
the summed ADC signal of the left HRS gas Cerenkov detec- 
tor, without a cut on the lead glass signal and after applying 
such lead glass electron and pion cuts. The spectrum from 
the right HRS is similar. Fig. [8] shows the distribution of the 
energy deposit in the two layers of the right HRS lead glass 
counters, without a Cerenkov cut, and after Cerenkov electron 
and pion cuts. 

Detailed PID analysis was done both before and during 
the experiment. The PID performance of each detector is 
characterized by the electron detection efficiency ij e and the 
pion rejection factor rj nxe i> defined as the number of pions 
needed to cause one pion contamination event. In the HRS 
central momentum range of 0.8 < po < 2.0 (GeV/c), the PID 
efficiencies for the left HRS were found to be 

o Gas Cerenkov: -q^^ > 770 at r\ e = 99.9%; 

o Lead glass counters: 77^,^ ~ 38 at rj e = 98%; 

o Combined: ?77r,rcj > 3 x 10 4 at? ?e = 98%. 

and for the right HRS were 

o Gas Cerenkov: r) n Ie \ = 900 at r} e = 99%; 

o Lead glass counters: ?/ x ,roj ~ 182 at r/ e = 98%; 

o Combined: rj^^j > 1.6 x 10 5 at ?/ e = 97%. 
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FIG. 7: (Color online) Summed ADC signal of the left HRS gas 
Cerenkov detector: without cuts, after lead glass counters electron 
cut and after pion cut. The vertical line shows a cut Yl ADCi > 400 
applied to select electrons. 
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FIG. 8: (Color online) Energy deposited in the first layer (preshower) 
vs that in the second layer (shower) of lead glass counters in the 
right HRS. The two blobs correspond to the spectrum with a tight 
gas Cerenkov ADC electron cut and with a pion cut applied. The 
lines show the boundary of the two-dimensional cut used to select 
electrons in the data analysis. 



analysis, and to calibrate experiment-specific detector con- 
stants. It also provided the possibility to apply conditions on 
the incoming data. The information from scaler events was 
used to extract beam charge and DAQ deadtime corrections. 



IV. THE POLARIZED TARGET 

Polarized 3 He targets are widely used at SLAC, DESY, 
MAINZ, MIT-Bates and JLab to study the electromagnetic 
structure and the spin structure of the neutron. There exist 
two major methods to polarize 3 He nuclei. The first one uses 
the metastable-exchange optical pumping technique l70ll . The 
second method is based on optical pumping 1 7 1 ] and spin ex- 
change If72ll . It has been used at JLab since 1998 If73ll . and 
was used here. 

The 3 He target at JLab Hall A uses the same design as the 
SLAC 3 He target The first step to polarize 3 He nuclei 
is to polarize an alkali metal vapor (rubidium was used at 
JLab as well as at SLAC) by optical pumping ifTlll with cir- 
cularly polarized laser light. Depending on the photon helic- 
ity, the electrons in the Rb atoms will accumulate at either 
the F = 3, m p = 3 or the F = 3, nip = —3 level (here 
F is the atom's total spin and mp is its projection along the 
magnetic field axis). The polarization is then transfered to the 
3 He nuclei through the spin exchange mechanism 17211 during 
collisions between Rb atoms and the 3 He nuclei. Under oper- 
ating conditions the 3 He density is about 10 20 nuclei/cm 3 and 
the Rb density is about 10 14 atoms/cm 3 . 

To minimize depolarization effects caused by the unpolar- 
ized light emitted from decay of the excited electrons, N2 
buffer gas was added to provide a channel for the excited 
electrons to decay to the ground state without emitting pho- 
tons l7lll . In the presence of N2, electrons decay through 
collisions between the Rb atoms and N2 molecules, which is 
usually referred to as non-radiative quenching. The number 
density of N2 was about 1% of that of 3 He. 



A. Target Cells 

The target cells used for this experiment were 25-cm lon£ 
pressurized glass cells with ~ 130-/J,m thick end windows. 



K. Data Acquisition System 

We used the CEBAF Online Data Acquisition (CODA) sys- 
tem l68ll for this experiment. In the raw data file, data from 
the detectors, the beamline equipment, and from the slow con- 
trol software were recorded. The total volume of data ac- 
cumulated during the two-month running period was about 
0.6 TBytes. Data from the detectors were processed using 
an analysis package called Experiment Scanning Program for 
hall A Collaboration Experiments (ESPACE) @. ESPACE 
was used to filter raw data, to make histograms for recon- 
structed variables, to export variables into ntuples for further 



TABLE III: Target cell characteristics. Symbols are: V p pumping 
chamber volume in cm 3 ; Vt target chamber volume in cm 3 ; Vtr 
transfer tube volume in cm 3 ; Vb total volume in cm 3 ; L tr transfer 
tube length in cm; no: 3 He density in amg at room temperature (1 
amg = 2.69 x 10~ 19 /cm 3 which corresponds to the gas density at 
the standard pressure and T = 0°C); lifetime is in hours. 
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Cell #2 
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3.9 


173.5 
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8.28 


44 


uncertainty 


1.5 
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0.020 
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The cell consisted of two chambers, a spherical upper cham- 
ber which holds the Rb vapor and in which the optical pump- 
ing occurs, and a long cylindrical chamber where the electron 
beam passes through and interacts with the polarized 3 He nu- 
clei. Two cells were used for this experiment. Figure [9] is a 
picture of the first cell with dimensions shown in mm. Ta- 
ble|^gives the cell volumes and densities. 



B. Target Setup 

Figure^|is a schematic diagram of the target setup. There 
were two pairs of Helmholtz coils to provide a 25 G main 
holding field, with one pair oriented perpendicular and the 
other parallel to the beamline (only the perpendicular pair is 
shown). The holding field could be aligned in any horizon- 
tal direction with respect to the incident electron beam. The 
coils were excited by two power supplies in the constant volt- 
age mode. The coil currents were continuously measured and 
recorded by the slow control system. The cell was held at 
the center of the Helmholtz coils with its pumping chamber 
mounted inside an oven heated to 170°C in order to vapor- 
ize the Rb. The lasers used to polarize the Rb were three 
30 W diode lasers tuned to a wavelength of 795 nm. The tar- 
get polarization was measured by two independent methods - 
the NMR (Nuclear Magnetic Resonance) 1 6 dl73l.| 75 jl and the 
EPR (Electro Paramagnetic Resonance) I58L I6CL VTjL 17611 po- 
larimetry. The NMR system consisted of one pair of pick-up 
coils (one on each side of the cell target chamber), one pair 
of RF coils and the associated electronics. The RF coils were 
placed at the top and the bottom of the scattering chamber, ori- 
ented in the horizontal plane, as shown in Fig. Ell The EPR 
system shared the RF coils with the NMR system. It consisted 
of one additional RF coil to induce light signal emission from 
the pumping chamber, a photodiode and the related optics to 
collect the light, and associated electronics for signal process- 
ing. 



FIG. 10: (Color online) Target setup overview (schematic). 



C. Laser System 



The laser system used during this experiment consisted of 
seven diode lasers - three for longitudinal pumping, three for 
transverse pumping and one spare. To protect the diode lasers 
from radiation damage from the electron beam, as well as 
to minimize the safety issues related to the laser hazard, the 
diode lasers and the associated optics system were located in 
a concrete laser hut located on the right side of the beamline at 
90°, as shown in Fig.|4] The laser optics had seven individual 
lines, each associated with one diode laser. All seven optical 
lines were identical and were placed one on top of the other 
on an optics table inside the laser hut. Each optical line con- 
sisted of one focusing lens to correct the angular divergence 
of the laser beam, one beam-splitter to linearly polarize the 
lasers, two mirrors to direct them, three quarter waveplates 
to convert linear polarization to circular polarization, and two 
half waveplates to reverse the laser helicity. Figure [TT1 shows 
a schematic diagram of one optics line. 

Under the operating conditions for either longitudinal or 
transverse pumping, the original beam of each diode laser was 
divided into two by the beam-splitter. Therefore there were a 
total of six polarized laser beams entering the target. The di- 
ameter of each beam was about 5 cm which approximately 
matched the size of the pumping chamber. The target was 
about 5 m away from the optical table. For the pumping of 
the transversely polarized target, all these laser beams went 
directly towards the pumping chamber of the cell through a 
window on the side of the target scattering chamber enclo- 
sure. For longitudinal pumping, they were guided towards the 
top of the scattering chamber, then were reflected twice and 
finally reached the cell pumping chamber. 
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FIG. 11: (Color online) Laser polarizing optics setup (schematic) for 
the Hall A polarized 3 He target. 



D. NMR Polarimetry 

The polarization of the 3 He was determined by measuring 
the 3 He Nuclear Magnetic Resonance (NMR) signal. The 
principle of NMR polarimetry is the spin reversal of 3 He nu- 
clei using the Adiabatic Fast Passage (AFP) ll77ll technique. 
At resonance this spin reversal will induce an electromagnetic 
field and a signal in the pick-up coil pair. The signal mag- 
nitude is proportional to the polarization of the 3 He and can 
be calibrated by performing the same measurement on a wa- 
ter sample, which measures the known thermal polarization 
of protons in water. The systematic error of the NMR mea- 
surement was about 3%, dominated by the error in the water 
calibration Il75ll. 



E. EPR Polarimetry 



quency shift <Wepr is proportional to the 3 He polarization 
in the pumping chamber. The 3 He polarization in the target 
chamber is calculated using a model which describes the po- 
larization diffusion from the pumping chamber to target cham- 
ber. The value of the EPR resonance frequency ^epr can also 
be used to calculate the magnetic field magnitude. The sys- 
tematic error of the EPR measurement was about 3%, which 
came mainly from uncertainties in the cell density and tem- 
perature, and from the diffusion model. 



E Target Performance 

The target polarizations measured during this experiment 
are shown in Fig. [2] Results from the two polarimetries are 
in good agreement and the average target polarization in beam 
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FIG. 12: Target polarization, starting June 1 of 2001, as measured by 
EPR and NMR polarimetries. 

was (40.0 ± 2.4)%. In a few cases the polarization measure- 
ment itself caused an abrupt loss in the polarization. This phe- 
nomenon may be the so-called "masing effect" l74ll due to 
non-linear couplings between the 3 He spin rotation and con- 
ducting components inside the scattering chamber, e.g., the 
NMR pick-up coils, and the "Rb-ring" formed by the rubid- 
ium condensed inside the cell at the joint of the two cham- 
bers. This masing effect was later suppressed by adding coils 
to produce an additional field gradient. 



In the presence of a magnetic field, the Zeeman splitting 
of Rb, characterized by the Electron-Paramagnetic Resonance 
frequency ^epr, is proportional to the field magnitude. When 
3 He nuclei are polarized (P « 40%), their spins generate a 
small magnetic field £?3 Ho of the order of « 0.1 Gauss, super- 
imposed on the main holding field Bh = 25 Gauss. During 
an EPR measurement ll76ll the spin of the 3 He is flipped by 
AFP, hence the direction of i?3 He is reversed and the change 
in the total field magnitude causes a shift in ^epr- This fre- 



V. DATA ANALYSIS 

In this section we present the analysis procedure leading 
to the final results in Section Ivfl We start with the analysis 
of elastic scattering, the A(1232) transverse asymmetry, and 
the check for false asymmetry. Next, the DIS analysis and 
radiative corrections are presented. Finally we describe nu- 
clear corrections which were used to extract neutron structure 
functions from the 3 He data. 
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A. Analysis Procedure 



spectrometer acceptance based on a transport model for the 



The procedure to extract the electron asymmetries from our 
data is outlined in Fig. ^] From the raw data one first ob- 
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FIG. 13: Procedure for asymmetry analysis. 

tains the helicity-dependent electron yield N ± using accep- 
tance and PID cuts. The efficiencies associated with these 
cuts are not helicity-dependent, hence are not corrected for in 
the asymmetry analysis. The yield is then corrected for the 
helicity-dependent integrated beam charge Q ± and the live- 
time of the DAQ system LT*. The asymmetry of the cor- 
rected yield is the raw asymmetry A raw . Next, to go from 
Araw to the physics asymmetries Aii and A±, four factors 
need to be taken into account: the beam polarization Pj,, the 
target polarization P t , the nitrogen dilution factor /n 2 due to 
the unpolarized nitrogen nuclei mixed with the polarized 3 He 
gas, and a sign based on the knowledge of the absolute state 
of the electron helicity and the target spin direction: 

A ^ = ± l^Kp (17) 

The results of the beam and the target polarization measure- 
ments have been presented in previous sections. The nitrogen 
dilution factor is obtained from data taken with a reference 
cell filled with nitrogen. The sign of the asymmetry is 
described by "the sign convention". The sign convention for 
parallel asymmetries was obtained from the elastic scattering 
asymmetry and that for perpendicular asymmetries was from 
the A(1232) asymmetry analysis, as will be described in 
Sections fV Bl and ly CI The physics asymmetries An and A±, 
after corrections for radiative effects, were used to calculate 
A\ and A2 and the structure function ratios g\jF\ and 32/^1 
using Eq. (IA22I — lA25l . Then the last step is to apply nuclear 
corrections in order to extract the neutron asymmetries and 
the structure function ratios from the 3 He results, as will be 
described in Section|VF| 



Although the main goal of this experiment was to provide 
precise data on the asymmetries, cross sections were also ex- 
tracted from the data. The procedure for the cross section 
analysis is outlined in Fig. ^] One first determines the ab- 
solute yield of e — 3 Hc inclusive scattering from the raw data. 
Unlike the asymmetry analysis, corrections need to be made 
for the detector and PID efficiencies and the spectrometer ac- 
ceptance. A Monte-Carlo simulation is used to calculate the 
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FIG. 14: Procedure for cross section analysis. 

HRS IrSOll with radiative effects taken into account. One then 
subtracts the yield of e — N scattering caused by the N2 nuclei 
in the target. The clean e — 3 Hc yield is then corrected for 
the helicity-averaged beam charge and the DAQ livetime to 
give cross section results. Using world fits for the unpolarized 
structure functions (form factors) of 3 He, one can calculate 
the expected DIS (elastic) cross section from the Monte-Carlo 
simulation and compare to the data. 



B. Elastic Analysis 

Data for e — 3 He elastic scattering were taken on a lon- 
gitudinally polarized target with a beam energy of 1.2 GeV. 
The scattered electrons were detected at an angle of 20°. The 
formalism for the cross sections and asymmetries are summa- 
rized in AnnendixlBl Results for the elastic asymmetry were 
used to check the product of beam and target polarizations, 
as well as to determine the sign convention for different 
beam-helicity states and target spin directions. 

The raw asymmetry was extracted from the data by 
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(18) 



Q+LT+ t Q-LT- 



with N^, and LT ± the helicity-dependent yield, beam 
charge and livetime correction, respectively. The elastic 
asymmetry is 



± 



fN 2 fQEPbPt 



(19) 



with / N , = 0.975 ± 0.003 the N 2 dilution factor determined 
from data taken with a reference cell filled with nitrogen, and 
Pb and Pt the beam and target polarization, respectively. A 
cut in the invariant mass \ W — A^3 Ho | < 6 (MeV) was used to 
select elastic events. Within this cut there are a small amount 
of quasi-elastic events and Jqe > 0.99 is the quasi-elastic 
dilution factor used to correct for this effect. 

The sign on the right hand side of Eq. H9i depends on the 
configuration of the beam half-wave plate, the spin precession 
of electrons in the accelerator, and the target spin direction. It 
was determined by comparing the sign of the measured raw 
asymmetries with the calculated elastic asymmetry. We found 
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FIG. 15: (Color online) Elastic parallel asymmetry results for the two 
HRS. The kinematics are E = 1.2 GeV and 9 = 20°. A cut in the 
invariant mass \W — M3 He | < 6 (MeV) was used to select elastic 
events. Data from runs with beam half-wave plate inserted are shown 
as triangles. The error bars shown are total errors including a 4.5% 
systematic uncertainty, which is dominated by the error of the beam 
and target polarizations. The combined asymmetry and its total error 
from w 20 runs are shown by the horizontal solid and dashed lines, 
respectively, as well as the solid circle as labeled 15^1 . 

that for this experiment the electron helicity was aligned to the 
beam direction during H+ pulses when the beam half-wave 
plate was not inserted. Since the electron spin precession in 
the accelerator can be well calculated using quantum electro- 
dynamics and the results showed that the beam helicity dur- 
ing H+ pulses was the same for the two beam energies used 
for elastic and DIS measurements, the above convention also 
applies to the DIS data analysis. 

A Monte-Carlo simulation was performed which took into 
account the spectrometer acceptance, the effect of the quasi- 
elastic scattering background and radiative effects. Results 
for the elastic asymmetry and the cross section are shown in 
Fig. and ^] respectively, along with the expected values 
from the simulation. The data show good agreement with the 
simulation within the uncertainties. 



C. A(1232) Transverse Asymmetry 

Data on the A(1232) resonance were taken on a trans- 
versely polarized target using a beam energy of 1.2 GeV. The 
scattered electrons were detected at an angle of 20° and the 
central momentum of the spectrometers was set to 0.8 GeV/c. 
The transverse asymmetry defined by Eq. dA15i was extracted 
from the raw asymmetry using Eq. Mil . A cut in the invariant 
mass \W - 1232| < 20 (MeV) was used to select A(1232) 
events. The sign on the right hand side of Eq. Ml\ depends on 
the beam half-wave plate status, the spin precession of elec- 
trons in the accelerator, the target spin direction, and in which 



FIG. 16: (Color online) Elastic cross section results for the two HRS. 
The kinematics were E = 1.2 GeV and 9 = 20°. A systematic er- 
ror of 6.7% was assigned to each data point, which was dominated 
by the uncertainty in the target density and the HRS transport func- 
tions |5S|. 
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FIG. 17: (Color online) Measured raw A(1232) transverse asymme- 
try, with beam half-wave plate inserted and target spin pointing to the 
left side of the beamline. The kinematics are E — 1.2 GeV, 9 = 20° 
and E' = 0.8 GeV/c. The dashed lines show the expected value 
obtained from previous 3 He data extrapolated in Q 2 . 



(left or right) HRS the asymmetry is measured. Since data 
from a previous experiment 117 311 in a similar kinematic region 
showed that Af < and > jH, Af can be used 
to determine the sign convention of the measured transverse 
asymmetries. The raw A(1232) transverse asymmetry mea- 
sured during this experiment was positive on the left HRS, 
as shown in Fig. ^] with the beam half-wave plate inserted 
and the target spin pointing to the left side of the beamline. 
Also shown is the expected value obtained from previous 3 He 
data extrapolated in Q 2 . Similar to the longitudinal configu- 
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ration, this convention applied to both the A(1232) and DIS 
measurements. 



D. False Asymmetry and Background 

False asymmetries were checked by measuring the asym- 
metries from a polarized beam scattering off an unpolarized 
12 C target. The results show that the false asymmetry was 
less than 2 x 10 -3 , which was negligible compared to the sta- 
tistical uncertainties of the measured 3 He asymmetries. To es- 
timate the background from pair production 7 — > e~ + e + , the 
positron yield was measured at x = 0.33, which is expected 
to have the highest pair production background. The positron 
cross section was found to be ss 3% of the total cross section 
at x = 0.33, and the positron contribution at x = 0.48 and 
x = 0.61 should be even smaller. The effect of pair produc- 
tion asymmetry is negligible compared to the statistical uncer- 
tainties of the measured 3 He asymmetries and is not corrected 
for in this analysis. 



TABLE IV: Total radiation length Xo and thickness d of the ma- 
terial traversed by incident (before interaction) and scattered (after 
interaction) electrons. The cell is made of glass GE180 which has 
Xo = 7.04 cm and density p — 2.77 g/cm 3 . The radiation length 
and thickness after interaction are given by left/right depending on 
by which HRS the electrons were detected. 



X 


0.33, 0.48 


0.61 


0.61 


e 


35° 


45° 


45° 


Cell 


#2 


#2 


#1 


Cell window (/im) 


144 


144 


132 


Xo (before) 


0.00773 


0.00773 


0.00758 


d (g/cm 2 , before) 


0.23479 


0.23479 


0.23317 


Cell wall (mm) 


1.44/1.33 


1.44/1.33 


1.34/1.43 


X (after) 


0.0444/0.0416 


0.0376/0.0354 


0.0356/0.0374 


d (g/cm 2 , after) 


0.9044/0.8506 


0.7727/0.7293 


0.7336/0.7687 



TABLE V: Internal radiative corrections to Aii Hc and A 



E. DIS Analysis 

The longitudinal and transverse asymmetries defined by 
Eq. (IA13i and ( IA15i for DIS were extracted from the raw 
asymmetries as 



X 


AA^ 3He (xlO" 3 ) 


AA^ 3h ° (x10~ 3 ) 


0.33 


-5.77 ± 0.47 


2.66 ± 0.03 


0.48 


-3.28 ±0.13 


1.47 ± 0.05 


0.61 


-2.66 ±0.15 


1.28 ± 0.07 



.4 



±- 



~fa 3 P b Pt 



(20) 



where the sign on the right hand side was determined by the 
procedure described in Sections fV Bl and lV CI The N2 dilution 
factor, extracted from runs where a reference cell was filled 
with pure N 2 , was found to be / Na 0.938 ± 0.007 for all three 
DIS kinematics. 



Radiative corrections were performed for the 3 He asym- 
metries ^4j| Hc and A^- c . We denote by A ohs the observed 

asymmetry, A Born the non-radiated (Born) asymmetry, AA lr 
the correction due to internal radiation effects and A^4 er the 
one due to external radiation effects. One has A Born = 
A obs + AA lr + AA er for a specific target spin orientation. 

Internal corrections were calculated using an improved ver- 
sion of POLRAD 2.0 17^1 . External corrections were calcu- 
lated with a Monte-Carlo simulation based on the procedure 
first described by Mo and Tsai jsoil. Since the theory of radia- 
tive corrections is well established l80ll . the accuracy of the 
radiative correction depends mainly on the structure functions 
used in the procedure. To estimate the uncertainty of both 
corrections, five different fits 18 ll l82t l83l l84l I85J1 were used 
for the unpolarized structure function F2 and two fits l86H87ll 
were used for the ratio R. For the polarized structure function 
gu in addition to those used in POLRAD 2.0 HHH], we fit to 
world g\ j Ff and g™ j F[ l data including the new results from 
this experiment. Both fits will be presented in Section lVIBI 
For 52 we used both g^ w and an assumption that g^ = 0. 
The variation in the radiative corrections using the fits listed 
above was taken as the full uncertainty of the corrections. For 



TABLE VI: External radiative corrections to ^4|| He and A x He . Errors 
are from uncertainties in the structure functions and in the cell wall 
thickness. 



X 


A^ r ' 3Hc (xlO- 3 ) 


AAl r ' 3Ho (xl0- 3 ) 


0.33 


-0.67 ±0.10 


-0.05 ±0.11 


0.48 


-1.16 ±0.15 


0.80 ± 0.46 


0.61 


-0.39 ± 0.03 


0.29 ± 0.04 



external corrections the uncertainty also includes the contribu- 
tion from the uncertainty in the target cell wall thickness. The 
total radiation length and thickness of the material traversed 
by the scattered electrons are given in Table llVl for each kine- 
matic setting. Results for the internal and external radiative 
corrections are given in Table IVland l VII respectively. 

By measuring DIS unpolarized cross sections and using the 
asymmetry results, one can calculate the polarized cross sec- 
tions and extract gi and 52 from Eq. dA5l > and iA6\ . We used 
a Monte-Carlo simulation to calculate the expected DIS un- 
polarized cross sections within the spectrometer acceptance. 
This simulation included internal and external radiative cor- 
rections. The structure functions used in the simulation were 
from the latest DIS world fits l83H87ll with the nuclear effects 
corrected l90ll . The radiative corrections from the elastic and 
quasi-elastic processes were calculated in the peaking approx- 
imati on 19111 using the world proton and neutron form factor 
data I92L l93L 19411 . The DIS cross section results agree with 
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the simulation at a level of 10%. Since this is not a dedicated 
cross section experiment, we obtained the values for gi and g 2 
by multiplying our g\jF\ and 172 / Fi results by the world fits 
for unpolarized structure functions F\ I83ll87ll . instead of the 
Fx from this analysis. 



F. From He to Neutron 



Eq. ( I22> was also used for extracting A 2 , <7™ / 'F™ and g 2 / F" 
from our 3 He data. The uncertainty in A" due to the uncer- 
tainties in F^^, in the correction for EMC effects, in A\ data 
and in P n „ is given in Table|X] Compared to the convolution 
approach I98ll used by previous 3 He experiments ll50ll5lll52il . 
in which only the first two terms on the right hand side of 
Eq. d2 1 1 are present, the values of A" extracted from Eq. (1221 
are larger by (1 — 2)% in the region 0.2 < x < 0.7. 



Properties of protons and neutrons embedded in nuclei are 
expected to be different from those in free space because of 
a variety of nuclear effects, including that from spin depo- 
larization, binding and Fermi motion, the off-shell nature of 
the nucleons, the presence of non-nucleonic degrees of free- 
dom, and nuclear shadowing and antishadowing. A coherent 
and complete picture of all these effects for the 3 He structure 
function g^ in the range of 10~ 4 < x ^ 0.8 was presented 
in l97ll . It gives 



5i Hc - P n g?+2P p £-0.0U[fi(x)-4g?(x)] 
+a{x)gUx)+b{x)gl{x) 



(21) 



where P n (P p ) is the effective polarization of the neutron 
(proton) inside 3 He l57ll . Functions a(x) and b(x) are Q 2 - 
dependent and represent the nuclear shadowing and antishad- 
owing effects. 

From Ea. JA12> . the asymmetry A\ is approximately the ra- 
tio of the spin structure function g\ and F\ . Noting that shad- 
owing and antishadowing are not present in the large x region, 
using Eq. <ED 

one obtains 



G. Resonance Contributions 



Since there are a few nucleon resonances with masses above 
2 GeV and our measurement at the highest x point has an in- 
variant mass close to 2 GeV, the effect of possible contribu- 
tions from baryon resonances were evaluated. This was done 
by comparing the resonance contribution to g" with that to 
F". For our kinematics at x = 0.6, data on the unpolarized 
structure function F2 and R l95ll show that the resonance con- 
tribution to F\ is less than 5%. The resonance asymmetry was 
estimated using the MAID model |9(| and was found to be 
approximately 0.10 at W = 1.7 (GeV). Since the resonance 
structure is more evident at smaller W, we took this value as 
an upper limit of the contribution at W = 2 (GeV). The res- 
onance contribution to our A™ and gf /F™ results at x = 0.6 
were then estimated to be at most 0.008, which is negligible 
compared to their statistical errors. 



A™ = 



F^[A\^-2^P P A\{1 



0.014 n 

2P D > 



The two terms 



0.056 



and 



PnF?{\ 



0.014 



0.056 ^ 



(22) 



— represent the corrections to A™ 

associated with the A(1232) component in the 3 He wave- 
function. Both terms cause A" to increase in the x range 
of this experiment, and to turn positive at lower values of x 
compared to the situation when the effect of the A(1232) is 
ignored. For F 2 and F 2 Hc , we used the world proton and 
deuteron F 2 data and took into account the EMC effects l90ll . 
We used the world proton asymmetry data for A\. The ef- 
fective nucleon polarizations P„ iP can be calculated using 
3 He wavefunctions constructed from N-N interactions, and 
their uncertainties were estimated using various nuclear mod- 
els HISS, giving 



P„ 



0.86 



+0.036 
0.02 



and P„ 



-0.028 



+0.009 
0.004 



(23) 



VI. RESULTS 



A. 3 He Results 



Results of the electron asymmetries for e — 3 He scattering, 
A|| Hc and A^ c , the virtual photon asymmetries A^ and 

A 2 Hc , structure function ratios .9i Hc /-Fi Hc and <?2 H °/^i He 
and polarized structure functions g'^ and g 2 Hc are given in 
Table IVIII Results for g^ were obtained by multiplying 
the gi^/F^ results by the unpolarized structure function 

F^ , which were calculated using the latest world fits of 
DIS data |H |H and with nuclear effects corrected l90ll . 
Results fo r and gA^ are show n in Fig. ^] along with 
SLAC HHEM and HERMES OH data. 
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FIG . 18: Results for the 3 He asymmetry A 1 Uc and the structure functions <7i He as a function of x, along with previous data from SLAC l5ll 



100] and HERMES IIOlll . Error bars of the results from this work include both statistical and systematic uncertainties. 



TABLE VII: Results for 3 He asymmetries Al Hc and A 3 2 Hc , structure function ratios gl Be /F^ 
functions <7! Hc and ;?2 Ho - Errors are given as ± statistical ± systematic. 



and 32 Ho /^i Hc > an d polarized structure 



(x) 


0.33 


0.47 


0.60 


(Q 2 ) (GeV/c) 2 


2.71 


3.52 


4.83 


JJ 


-0.020 ± 0.005 ±0.001 


-0.012 ±0.005 ±0.000 


0.007 ± 0.007 ± 0.001 


a bc 


0.000 ± 0.010 ± 0.000 


0.016 ±0.008 ±0.001 


-0.010 ±0.016 ±0.001 




-0.024 ± 0.006 ±0.001 


-0.019 ±0.006 ±0.001 


0.010 ±0.009 ±0.001 


A 3 ^ 


-0.004 ±0.014 ±0.001 


0.020 ±0.012 ±0.001 


-0.013 ±0.023 ±0.001 


si H 7^i He 


-0.022 ±0.005 ±0.001 


-0.008 ± 0.008 ± 0.001 


0.003 ± 0.009 ± 0.001 


„ 3 He/p 3 He 

92 1 f\ 


0.010 ± 0.036 ± 0.002 


0.050 ± 0.022 ± 0.003 


-0.028 ± 0.038 ± 0.002 




-0.024 ± 0.006 ±0.001 


-0.004 ± 0.004 ± 0.000 


0.001 ± 0.002 ± 0.000 


9 > 


0.011 ±0.039 ±0.001 


0.026 ± 0.012 ± 0.002 


-0.006 ± 0.009 ± 0.001 



B. Neutron Results 



Results for the neutron asymmetries A™ and , structure function ratios /F" and g% /F™ and polarized structure functions 
<7™ and are given in Table I Villi 



TABLE VIII: Results for the asymmetries and spin structure functions for the neutron. Errors are given as ±statistical ±systematic. 



(x) 


0.33 


0.47 


0.60 


(Q 2 ) (GeV/c) 2 


2.71 


3.52 


4.83 


A" 


-0.048 ± 0.024±g;gl| 


-0.006 ± 0.027±g:gif 


0.175 ± 0.048±g:g|| 


Al 


-0.004 ± 0.063±£ 


0.117 ± 0.055±g;gaf 


-0.034 ±0.1241°;°^ 




-0.043 ± 0.022±g;gg| 


0.040 ± 0.035±g;gii 


0.124 ±0.045to.ol? 


gS/F? 


0.034 ± 0.153±g;g^g 


0.207 ±0.103tg;^ 2 


-0.190 ±0.204l° ° 2 ^ 


gi 


-0.012 ± 0.006lo;o°3 


0.005 ± D.004±g:ggl 


0.006 ± 0.002To.ooi 


gZ 


0.009 ± 0.043±o;go3 


0.026 ± 0.013to.So3 


-0.009 ± 0.009l°;gol 
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The A", g" /F" and g" results are shown in Fig. [T9ll20l and 
ETI respectively. In the region of x > 0.4, our results have 
improved the world data precision by about an order of mag- 
nitude, and will provide valuable inputs to parton distribution 
function (PDF) parameterizations. Our data at x = 0.33 are in 
good agreement with previous world data. For the A™ results, 
this is the first time that the data show a clear trend that A™ 
turns to positive values at large x. As x increases, the agree- 
ment between the data and the predictions from the constituent 
quark models (CQM) becomes better. This is within the ex- 
pectation since the CQM is more likely to work in the valence 
quark region. It also indicates that A" will go to higher values 
at x > 0.6. However, the trend of the A™ results does not 
agree with the BBS and LSS(BBS) parameterizations, which 
are from leading-order pQCD analyses based on hadron he- 
licity conservation (HHC). This indicates that there might be 
problem in the assumption that quarks have zero orbital angu- 
lar momentum which is used by HHC. 
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FIG. 19: Our A" results along with theoretical predictions and 
previous world data obtained from polarized 3 He targets 15(1 l5ll 
1531 . Curves: predictions of A™ from SU(6) symmetry (x axis 
at zero) Q| , constituent quark model (shaded band) l2al . statisti- 
cal model at Q 2 =4 (GeV/c) 2 (long-dashed) l42ll . quark-hadron 
duality using two different SU(6) breaking mechanisms (dash-dot- 
dotted and dash-dot-dot-dotted), and non-meson cloudy bag model 
(dash-dotted) EH; predictions of g"/F? from pQCD HHC based 
BBS parameterization at Q 2 =4 (GeV/c) 2 (higher solid) 13011 and 
LSS(BBS) parameterization at Q 2 =4 (GeV/c) 2 (dashed) ^J, 
LSS 2001 NLO polarized parton densities at Q 2 = 5 (GeV/c) 2 
(lower solid) ^ and chiral soliton models at Q 2 = 3 (GeV/c) 2 
(long dash-dotted) and @] at Q 2 = 4.8 (GeV/c) 2 (dotted). 

The sources for the experimental systematic uncertainties 
are listed in Table IIXI Systematic uncertainties for the A™ re- 
sults include that from experimental systematic errors, uncer- 
tainties in internal radiative corrections AA™' tr and external 
radiative corrections AA™' er as derived from the values in Ta- 



FIG. 20: Results for g"/F™ along with previous world data from 
SLAC j3 The curves are the prediction for gi/F" from the 
LSS 2001 NLO polarized parton densities at Q 2 = 5 (GeV/c) 2 EUl . 



the E155 experimental fit at Q 



5 (GeV/cr (long dash-dot- 



dotted) 1 53] and the new fit as described in the text (long dash-dot- 
dot-dotted). 
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FIG. 21: Results for along with previous world data from 
SLAC and HERMES 



TABLE IX: Experimental systematic errors for the A" result. 



source 


error 


Beam energy Eb 


AE b /E b < 5 x 10" 4 


HRS central momentum po 


AE e /E e < 5 x 10" 4 Q23J 


HRS central angle do 


A6» < 0.1° QQ4J 


Beam polarization Pf, 


AP b /P b < 3% 


Target polarization Pt 


APt/Pt < 4% 


Target spin direction at 


Aa t < 1° 
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bles IVland IVII and that from nuclear corrections as described 
in Section IVF1 Table 1X1 gives these systematic uncertainties 
for the A™ results along with their statistical uncertainties. 
The total uncertainties are dominated by the statistical uncer- 
tainties. 



TABLE X: Total uncertainties for A 



(x) 


0.33 


0.47 


0.60 


Statistics 


0.024 


0.027 


0.048 


Experimental syst. 


0.004 


0.003 


0.004 


AA »,ir 


0.012 


0.013 


0.015 




0.002 


0.002 


0.003 


Fl Fi 


0.006 


0.008 


+0.005 
-0.010 


EMC effect 


0.001 


0.000 


0.009 




0.001 


0.005 


0.011 


P P 


+0.005 


+0.009 


+0.018 


-0.012 


-0.020 


-0.037 



We used five functional forms, x a P n (x)(l + (3/Q 2 ), to fit 
our q" /Fi" results combined with data from previous experi- 
ments I25j.l53ll . Here P n is the n*' l -order polynomial, n = 1,2 
for a finite a or n = 1, 2, 3 if a is fixed to be 0. The total num- 
ber of parameters is limited to ^ 5. For the Q 2 -dependence of 
gi/Fj, we used a term 1 + /3/Q 2 as in the E155 experimental 
fit 15311 . No constraints were imposed on the fit concerning 
the behavior of g\/F\ as x — > 1. The function which gives 
the smallest \ 2 value is g\ / F[ L = (a + bx + cx 2 ) (1 + /3/Q 2 ). 
The new fit is shown in Fig.|20] Results for the fit parameters 
are given in Table IXII and the covariance error matrix is 

" 1.000 -0.737 0.148 0.960 " 

-0.737 1.000 -0.752 -0.581 

e ~ 0.148 -0.752 1.000 -0.039 ' 

0.960 -0.581 -0.039 1.000 



TABLE XII: Result of the fit gf/Ff = x a (a + bx)(l + (3/Q 2 ). 



a = 0.813 ±0.049 
a = 1.231 ±0.122 
6= -0.413 ±0.216 
(3 = 0.030 ±0.124 



Figures|22]and|23]show the results for AV^ and xgQ , respec- 
tively. The precision of our data is comparable to the data 
from E155x experiment at SLAC 110211 . which is so far the 
only experiment dedicated to measuring 52 with published re- 
sults. 

To evaluate the matrix eleme nt we combined our g% re- 
sults with the E155x data I102ll . The average Q 2 of the E155x 
data set is ab out 5 (GeV/c) 2 . Following a similar procedure as 
used in Ref. Il02il . we assumed that gi (x, Q 2 ) is independent 
of Q 2 and g 2 oc (1 - x) m with m = 2 or 3 for x > 0.78 be- 
yond the measured region of both experiments. We obtained 
from Eq. l|6) 

dl = 0.0062 ±0.0028. (24) 

Compared to the value published previously Il02t . the uncer- 
tainty on c?2 nas been improved by about a factor of two. The 
large decrease in uncertainty despite the small number of our 
data points arises from the x 2 weighting of the integral which 
emphasizes the large x kinematics. The uncertainties on the 
integrand has been improved in the region x > 0.4 due to our 
#2 results at the two higher x points being more precise than 
that of E155x. While a negative value was predicted by lattice 
QCD (HI and most other models OEQ, the new result 
for g?2 suggests that the higher twist contribution is positive. 



TABLE XI: Result of the fit g" / F" = (a ± bx + cx 2 )(l + f3/Q 2 



a = -0.049 ± 0.052 
b= -0.162 ±0.217 
c = 0.698 ± 0.345 
(3 = 0.751 ±2.174 



Similar fits were performed to the proton world 
data E1IHB and function g*/Ff = x a {a+bx){l+f}/Q 2 ) 
was found to give the smallest \ 2 value. The new fit is shown 
in Fig. [2] of Section Hi Gl Results for the fit parameters are 
given in Table IXlfl and the covariance error matrix is 
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FIG. 22: Results for A% along with the best previous world 
data Il02ll . The curve gives the twist-2 contribution at Q 2 = 
4 (GcV/c) 2 calculated using the E155 experimental fit 1 5 
ofEq.0- 
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FIG. 23: Results for xg% along with the best previous world 
data Il02ll . The curve gives the twist-2 contribution at Q 2 — 
4 (GeV/c) 2 calculated using the El 55 experimental fit IBUl and g^ w 
ofEq. 0. 



C. Flavor Decomposition using the Quark-Parton Model 

Assuming the strange quark distributions s(x), s(x), As(x) 
and As(x) to be small in the region x > 0.3, and ignoring 
any Q 2 -dependence of the ratio of structure functions, one 
can extract polarized quark distribution functions based on the 
quark-parton model as 



Ag\(A + R du ) g?(l + 4R 



15Ff 



15F? 



and 



Ad + Ad 4<?™(1 + AR du ) g%(4 + R 



d + d 



du 



15F[R 



da, 



(25) 



(26) 



with R du = (d +_d)/(u +_u). Results for (Am + A u)/(u+ 
u) and (Ad 4- Ad)/(d + d) are given in Table Ixml As in- 
puts we used our own results for g™/F™, the world data on 
g\/F^ l58ll . and the ratio R du extracted f rom proton and 
deuteron unpolarized structure function data Il05ll . In a sim- 
ilar manner as for Eq.( I25> and J26i and ignoring nuclear ef- 
fects, one can also add the world data on g-^/F-J 1 to the fit- 
ted data set and extract these polarized quark distributions. 
The results are, however, consistent with those given in Ta- 
ble ED an d nave ver y similar error bars because the data on 
the deuteron in general have poorer precision than the data on 
the proton and the neutron data from this experiment. The 
results presented here have changed compared to the values 
published previously in Ref. I15| due to an error discovered 
in our fitting of R du from Ref. Il05ll . The analysis procedure 
is consistent with what was used in Ref. fl5ll . 

Figure [24J shows our results along with semi-inclusive 
data on (Aq + Aq)/ (q + q) obtained from recent results 
for Aq and Aq llOfjll by the H ERM ES collaboration, and 
the CTEQ6M unpolarized PDF lH07il . To estimate the ef- 
fect of the s and s contributions, we used tw o unp olar- 
ized PDF sets, CTEQ6M (Tol and MRST2001 Il08ll . and 



TABLE XIII: Results for the polarized quark distributions. The three 
uncertainties are those due to the gi/Fi statistical error, gi/FJ 1 
systematic uncertainty and the uncertainties of the g\ / Ff data, the 



R du fit and the correction for s and c quark contributions. 


(x) 


(Au + Au)/(u + u) 


(Ad + Ad) /(d + d) 


0.33 
0.47 
0.60 


0.545 ± 0.004 ± 0.002±£;°25 
0.649 ± 0.006 ± 0.002±g;§|| 
0.728 ± 0.006 ± 0.002^°;^* 


-0.352 ± 0.035 ± 0.0141°;°^ 
-0.393 ± 0.063 ± 0.020to;oS 
-0.440 ± 0.092 ± 0.0351^-^2 
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FIG. 24: Results for (Am + Au)/(u + u) and (Ad + Ad)/(d + d) 
in the quar k-par ton model, compared with se mi-inc lusive data from 
HERMES flOr! and CTEQ unpolarized PDfJT07|) as described in 
the text, the RCQM predictions (dash-dotted) [28], predictions from 



LSS 2001 NLO polarized parton densities at Q 



5 (GeV/cr 



(solid) 14 111 , the statistical model at Q z = 4 (GeV/c) 2 (long- 
dashed) 1431. the pQCD-based predictions with the HHC constraint 
(dashed) 13111 . the duality model using two different SU(6) break- 
ing mechanisms (dash-dot-dotted and dash-dot-dot-dotted) 14711 . and 
predictions from chiral soliton model at Q 2 = 4.8 (GeV/c) 2 (dot- 
ted) l44ll . The error bars of our data include the uncertainties given 
in Table IXIIll The shaded band near the horizontal axis shows the 
difference between Aqv / qv and (Aq + Aq) /(q + q) that needs to 
be added to the data when comparing with the RCQM calculation. 



three polarized PDF sets, AAC2003 fiol . BB2002 fTToft and 
GRSV2000 ll 1 ill . Fo r c and c c ontributions we used the two 
unpolarized PDF sets fl07ifl03l and the positivity conditions 
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that |Ac/c| ^ 1 and |Ac/c| ^ 1. To compare with the 
RCQM predictions, which are given for valence quarks, the 
difference between Aqy / qv and ( Aq + Aq) / (a + q ) wa s es- 
timated using the two unpola rized PDF sets Il07l Il08ll and 
the three polarized PDF sets Il09l ll lOt ll 1 ill and is shown 
as the shaded band near the horizontal axis of Fig. EH Here 
qv{Aqv) is the unpolarized (polarized) valence quark distri- 
bution for u or d quark. Results shown in Fig. [24] agree well 
with the predictions from the RCQM I2H1 and the LSS 2001 
NLO polarized parton densities 14111 . The results agre e rea- 
sonably well with the statistical model calculation 14211 . But 
results for the d quark do not agree with the predictions from 
the leading-order pQCD LSS(BBS) parameterization 13111 as- 
suming hadron helicity conservation. 



or antiquark inside the target nucleon and transfers a large 
fraction of its energy and momentum via exchanged photons. 
In the single photon exchange approximation, the electron in- 
teracts with the target nucleon via only one photon, as shown 
in Fig. |25]||6|], and probes the quark structure of the nucleon 
with a spatial resolution determined by the four momentum 
transfer squared of the photon Q 2 = —q 2 . Moreover, if a po- 
larized electron beam and a polarized target are used, the spin 
structure of the nucleon becomes accessible. In the following 




q = (v, q) 



< 



VII. CONCLUSIONS 

We have presented precise data on the neutron spin 
asymmetry A™ and the structure function ratio g™ / F™ in the 
deep inelastic region at large x obtained from a polarized 
3 He target. These results will provide valuable inputs to the 
QCD parameterizations of parton densities. The new data 
show a clear trend that A™ becomes positive at large x. Our 
results for A\ agree with the LSS 2001 NLO QCD fit to the 
previous data and the trend of the x-dependence of A\ agrees 
with the hyperfine-perturbed RCQM predictions. Data on 
the transverse asymmetry and structure function A% and 
were also obtained with a precision comparable to the best 
previous world data in this kinematic region. Combined with 
previous world data, the matrix element was evaluated and 
the new value differs from zero by more than two standard 
deviations. This result suggests that the higher twist con- 
tribution is positive. Combined with the world proton data, 
the polarized quark distributions (Au + Au)/(u + u) and 
(Ad + Ad) /(d + d) were extracted based on the quark parton 
model. While results for (Au + Au) /(u + u) agree well with 
predictions from various models and fits to the previous data, 
results for (Ad + Ad) /(d + d) agree with the predictions 
from RCQM and from the LSS 2001 fit, but do not agree 
with leading order pQCD predictions that use hadron helicity 
conservation. Since hadron helicity conservation is based on 
the assumption that quarks have negligible orbital angular 
momentum, the new results suggest that the quark orbital 
angular momentum, or other effects beyond leading-order 
pQCD, may play an important role in this kinematic region. 



APPENDIX A: FORMALISM FOR ELECTRON DEEP 
INELASTIC SCATTERING 

The fundamental quark and gluon structure of strongly in- 
teracting matter is studied primarily through experiments that 
emphasize hard scattering from the quarks and gluons at suf- 
ficiently high energies. One important way of probing the dis- 
tribution of quarks and antiquarks inside the nucleon is elec- 
tron scattering, where an electron scatters from a single quark 



FIG. 25: (Color online) Electron scattering in the one-photon ex- 
change approximation. 



we denote the incident electron energy by E, the energy of the 
scattered electron by E' thus the energy transfer of the photon 
is v = E — E', and the three-momentum transfer from the 
electron to the target nucleus by q. 



1. Structure Functions 

In the case of unpolarized electrons scattering off an unpo- 
larized target, the differential cross-section for detecting the 
outgoing electron in a solid angle d£l and an energy range 
(E 1 , E' + dE') in the laboratory frame can be written as 

d 2 a /der\ 
dttdE' ~ \dnJMott ' 

-F 2 (x 1 Q 2 ) + ^-F 1 (x 1 Q 2 )t^ 29 - 



(Al) 



where 8 is the scattering angle of the electron in the laboratory 
frame. The four momentum transfer Q 2 is given by 



Q 2 = AEE' sin 2 - , 

and the Mott cross section, 

(da\ a 2 cos 2 § a 2 cos 2 | E' 

\dn) Matt ~ 4E 2 sin 4 § ~ 



E 



(A2) 



(A3) 



with a the fine structure constant, is the cross section for scat- 
tering relativistic electrons from a spin-0 point-like infinitely 
heavy target. Fi(x,Q 2 ) and F2(x,Q 2 ) are the unpolarized 
structure functions of the target, which are related to each 
other as 



Fi(x,Q 2 ) 



F 2 (x,Q 2 )(l+ 1 2 ) 
2x(l + R(x,Q 2 ) 



(A4) 



with 7 2 = (2Mx) 2 /Q 2 . Here R is defined as R = a L /a T 
with <tl and ot the longitudinal and transverse virtual photon 
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cross sections, which can also be expressed in terms of F\ and 
F 2 . 

Note that for a nuclear target, there exists an alternative per 
nucleoli definition (e.g. as used in Ref. |83|) which is 1/A 
times the definition used in this paper, here A is the number 
of nucleons inside the target nucleus. 

A review of doubly polarized DIS was given in Ref. fill . 
When the incident electrons are longitudinally polarized, the 
cross section difference between scattering off a target with its 
nuclear (or nucleon) spins aligned anti-parallel and parallel to 
the incident electron momentum is 

d 2 (T|4 d 2 CTj^. 4a 2 E' 



andE' dndE' vEQ 2 

x \{E + E' cos%i(x, Q 2 ) - 2Mxg 2 (x, Q 2 



.(AS) 



where gi(x,Q 2 ) and g2{x,Q 2 ) are the polarized structure 
functions. If the target nucleons are transversely polarized, 
then the cross section difference is given by 



d 2 <r T ^ d 2 <r T< 



'2 



andE' dridE' 



\u 2 E 
vEQ 2 

\2 



Slllt 



2E 

gi(x,Q z ) + —g 2 (x,( 



(A6) 



4. Virtual Photon-Nucleon Asymmetries 

Virtual photon asymmetries are defined in terms of a helic- 
ity de comp osition of the virtual photon absorption cross sec- 
tions ll 18tl . For the absorption of circularly polarized virtual 
photons with helicity ±1 by longitudinally polarized nucle- 
ons, the longitudinal asymmetry A\ is defined as 



A!(x,Q 2 ) 



cr l/2 — Q"3/2 
Cl/2 + C3/2 



(A7) 



where cr 1 / 2 (3/2) is the total virtual photo-absorption cross sec- 
tion for the nucleon with a projection of 1/2(3/2) for the total 
spin along the direction of photon momentum. 
A 2 is a virtual photon asymmetry given by 



A 2 {x,Q z ) 



2\ — 



2(7 LT 



f 1/2 + C3/2 



(A8) 



where <tlt describes the interference between transverse and 
longitudinal virtual photon-nucleon amplitudes. Because of 
the positivity limit, A 2 is usual ly sm all in the DIS region and 
it has an upper bound given by 



ly sma 

Eli 



A 2 {x,Q*)^\l?- l + A^Q 2 ) 



(A9) 



2. Bjorken Scaling and Its Violation 

A remarkable feature of the structure functions F\, F 2 
and g 2 is their scaling behavior. In the Bjorken limit 11 1 3ll 
(Q 2 — ► 00 and v — > 00 at a fixed v alue o f x), the structure 
functions become independent of Q 2 ill 14ll . Moreover, in this 
limit ctl vanishes |0], hence R = and Eq. dA4> redu ces t o 
£2(2:) = 2xFi(x), known as the Callan-Gross relation II 1511 . 

At finite Q 2 , the scaling of structure functions is violated 
due to the radiation of gluons by both initial and scattered 
quarks. These gluon radiative corrections cause a logarith- 
mic Q 2 -dependence to the struct ure f unctions, which has 
been verified by experimental data [1 16] and can be precisely 
calculated in pQCD using the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations Ill7ll . 



These two virtual photon asymmetries, depending in gen- 
eral on x and Q 2 , are related to the nucleon structure functions 

51(2, Q 2 ), g 2 (x,Q 2 ) andFx^Q 2 ) via 



Ai{x,Q 2 ) 



gi(x,Q 2 )-~/ 2 g 2 (x,Q 2 ) 



and 



M{x 7 Q 2 ) 



F 1 (x,Q 2 ) 



gi {x,Q 2 )+g 2 {x,Q 2 ) 



(A10) 



Fi(x,Q 2 ) 



At high Q 2 , one has j 2 -C 1 and 



Ai(x,Q 2 ) 



9i(x,Q 2 ) 
Fi(x,Q 2 ) 



(All) 



(A12) 



3. From Bjorken Limit to Finite Q 2 using the Operator 
Product Expansion 

In order to calculate observables at finite values of Q 2 , a 
method called the Operator Product Expansion (OPE) Q can 
be applied to DIS which can separate the non-perturbative 
part of an observable from its perturbative part. In the OPE, 
whether an operator is perturbative or not is characterized by 
the "twist" of the operator. At large Q 2 the leading twist t = 2 
term dominates, while at small Q 2 higher-twist operators need 
to be taken into account, which are sensitive to interactions 
beyond the quark-parton model, e.g., quark-gluon and quark- 
quark correlations 0] ■ 



In QCD the asymmetry A\ is expected to have less 
Q 2 -dependence than the structure functions themselves be- 
cause of the similar leading order Q 2 -evolution behavior of 
gi(x,Q 2 ) and Fi(x,Q 2 ). Existing data on the proton and 
the neutron asymmetries A\ and A™ indeed show little Q 2 - 
dependence 15311 . 



5. Electron Asymmetries 

In an inclusive experiment covering a large range of excita- 
tion energies the virtual photon momentum direction changes 
frequently, and it is usually more practical to align the target 
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spin longitudinally or transversely to the incident electron di- 
rection than to the momentum of the virtual photon. The vir- 
tual photon asymmetries can be related to the measured elec- 
tron asymmetries through polarization factors, kinematic vari- 
ables and the ratio R defined in Section lATI The longitudinal 
electron asymmetry is defined by Il2fil 



An 



gift - a n 
a n + a n 

(1 - e)M 3 
(1 - eR)vF l 



(E + E' cos %i 



Q 2 



-92 



(A13) 



where a^(a^) is the cross section of scattering off a longitu- 
dinally polarized target, with the incident electron spin aligned 
anti-parallel (parallel) to the target spin, and e is the magnitude 
of the virtual photon's longitudinal polarization: 



= l + 2(l + l/7 2 )tan 2 (6»/2) 



(A14) 



Similarly the transverse electron asymmetry is defined for a 
target polarized perpendicular to the beam direction as 



<7j^ - 

(1 - e)E'M 3 



(1 - eR)vFi 



2E 



9i 



-<I2 



cos0, (A15) 



where cri^(a^) is the cross section for scattering off a trans- 
versely polarized target with incident electron spin aligned 
anti-parallel (parallel) to the beam direction, and the scattered 
electrons being detected on the same side of the beam as that 
to which the target spin is pointing Jl5ll . The electron asym- 
metries can be written in terms of A\ and A2 as 



A,, = D(A 1 +r 1 A 2 ) 



(A16) 



and 



A± = d(A 2 -^A l ) , (A17) 
where the virtual photon polarization factor is given by 



D = 



1 + eR 



(A18) 



with y = v I ' E the fractional energy loss of the incident elec- 
tron. The remaining kinematic variables are given by 

V = (eVQ*)/(E-E'e) , (A19) 
£ = 77(1 + e)/(2e) and (A20) 
d = Dy/2e/(l + e) . (A21) 



6. Extracting Polarized Structure Functions from 
Asymmetries 

From Eq. dA16b and ( IA17I the virtual photon asymmetries 
A\ and A 2 can be extracted from measured electron asymme- 
tries as 



A x 



1 



D(l + vZ) 11 d(l + V 



A 



(A22) 



and 



A 2 = 



D(i + nO' 



-An 



d(l + tf) 



-A 1 



(A23) 



If the unpolarized structure functions F\ (x, Q 2 ) and R(x, Q 2 ) 
are known, then the polarized structure functions can b e ex- 
tracted from measured asymmetries Am and A± as 1 1 12] 



and 



9i(x,Q 2 ) 



92(x,( 



with D' given by 



Fi(x,Q 2 ) 



D' 



A\\ + Aj_taa(e/2) 



(A24) 



yF^x.Q 2 ) 
2D' sin (9 

, E + E' cos 6 
Ax ^ A\\ sm( 



D' 



(l-6)(2-y) 
y[l + eR{x,Q 2 ) 



,(A25) 



(A26) 



APPENDIX B: FORMALISM FOR E- 3 HE ELASTIC 
SCATTERING 

The cross section for electron elastic scattering off an un- 
polarized 3 He target can be written as 



da \ u (JMott f 2 



(—) 

Vdsi/ 



+ i w-[i|i- ta ° 2W2) l F ™ (<3) } <B1) 

where t = Q 2 /{AM 2 ) = u/(2M t ) is the recoil factor, M t 
is the target ( 3 He) mass, Q 2 is calculated from Eq. iA2\ . q is 
the three momentum transfer, /i3H e is the 3 He magnetic mo- 
ment, and F c and F m are the 3 He charge and magnetic form 
factors, which have been measured to a good precision 1 1 2. ill . 
The Mott cross section <JMott for a target of charge Z can be 
written as 



O Mott 



\ ail/ Mott 



Z 2 a 2 cos 2 (6>/2) E' 
4£ 2 sin 4 ((9/2) ~E 



with E' the energy of the outgoing electrons: 

E 



E' = 



I 4. IE. - 2 e 

1 + M T 2 



(B2) 



(B3) 



The elastic cross section for a polarized target can be writ- 
ten as 1H2I 



/ da \ h 



(da\ u 

\dn) 



hA(d*,0*,E,9,Q 2 ) , (B4) 



where h is the helicity of the incident electron beam, 
A(6* , </>*, E, 9, Q 2 ) describes the helicity-dependent cross 
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FIG. 26: (Color online) Polar and azimuthal angles of the target spin. 



and 



2J2t(1 + t)E' 
Rtv = - E — (ZF c )(fis Se F m ) , (Bll) 



The elastic asymmetry, defined by 



da \h=+l 



A 



el _ VdOdB' 



'_da_ 

, dQdE' 



'_da_\' l =+ 1 



V dfid-E' ) 



ft=T> ( B12 ) 



section, 0* is the polar angle and </>* is the azimuthal angle 
of the nucleon spin direction, as shown in Fig. [26] We write 
them expli citly for a target with spin parallel to the beam di- 
rection as Il22ll 



cos6>* = (E - E' cos 0)/\q\ 
cj)* = 0. 



(B5) 
(B6) 



The helicity-dependent part of the cross section can be written 

as 



— —&Mott 



(v t >Rt>(Q 2 ) cos 9* 



/da\ h =+ 1 /do-y^- 1 

+V TL >Rtl>(Q 2 ) sin I 
with kinematic factors 



Q 2 ,9 
55 tan 2V|F +tan 2 



and 



V TL > 



■ tan — 



(B8) 



(B9) 



V2\q\ 2 2 

Rt>, Rtl 1 can be related to the 3 He form factors F c , F m as: 
2tE' 



Rt> = — (M 3 Hcfm) 2 



(BIO) 



can therefore be calculated from Eq. JB 1 i and JB7> as 



Af=-(l-r) 

(Vt'Rt' cos 6* + Vtl'Rtl' sin 9* cos^* 

\ \q\ 2 c 2M 2 \2\q\ 2 LclU 1) r m] 



(B13) 
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